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Abstract 

We derive the exact expansion, to 0{rs), of the energy of the high-density spin-polarized two- 
dimensional uniform electron gas, where is the Seitz radius. 
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The three-dimensional uniform electron gas is a ubiquitous paradigm in solid-state physic^ 
and quantum chemistry^ and has been extensively used as a starting point in the development 
of exchange-correlation density functional in the framework of density- functional theoryPXhe 
two-dimensional version of the electron gas has also been the object of extensive researclJ^ 
because of its intimate connection to two-dimensional or quasi-two-dimensional materials, 
such as quantum dots.l^ 

The two-dimensional gas (or 2-jellium) is characterized by a density p = p^- + p^, where 
p^ and are the (uniform) densities of the spin-up and spin-down electrons, respectively. 
In order to guarantee its stability, the electrons are assumed to be embedded in a uniform 
background of positive chargePWe will use atomic units throughout. 

It is known from contributions by numerous workerd^^^ that the high-density {i.e. small-r^) 
expansion of the energy per electron (or reduced energy) in 2-jellium is 

^(^- C) = + + ^o(C) + ^.(C) Ts Inr, + 0(r,), (1) 

where = (vrp)^"*^^^ is the Seitz radius, and 

C = (2) 
P 

is the relative spin polarization.lSl Without loss of generality, we assume pi < pf, i.e. ( G [0, 1]. 

The first two terms of the expansion ([T]) are the kinetic and exchange energies, and their 
sum gives the Hartree-Fock (HF) energy. The paramagnetic {( = 0) coefficients are 

^-2(0) = +1, (3) 
^-i(0) = (4) 

and their spin-scaling functions are 

^-'^^^ = ^ = 2 ' 

In this Brief Report, we show that the next two terms, which dominate the expansion of the 
reduced correlation energy,'^^! can also be obtained in closed form for any value of the relative 
spin polarization (. 
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The logarithmic coefficient ^^(C) can be obtained by a Gell-Mann-Brueckner resummatiori^ 
of the most divergent terms in the infinite series in ([T]), and this yield^ 

1 r°° [ / 7, \ / 1, \ "1^ 



where 



and 



12y27r 
R(u) 



1 - 



(7) 



(8) 



(9) 



is the Fermi wave vector associated with the spin-up and spin-down electrons, respectively. 
After an unsuccessful attempt by Zia,'II'the paramagnetic {( = 0) and ferromagnetic {( = 1) 
values, 

£^(0) = ~V2 ~ ^) = -0-0863136 . . . , 
Ml) = ^.HO) = -i(^-l)= -0,0162582..., 
were found by Rajagopal and KimbalP^ and the spin-scaling function. 



T.(C) 



+ h + 3 



F{k^,k^) + F{k^,k^) 
10-37r 



(10) 
(11) 

(12) 



was obtained 30 years later by Chesi and Giuliani.'!^ The explicit expression for F[x, y) is 



F(x, y) = 4(x + y) — Tlx — AxE ( 1 ^ ) + 2x^K(a;, y) 



x" 



where 



ti{x,y) 



(x^ — y"^) ^1"^ arccos(?//x), x <y^ 



(13) 



(14) 



(y^ — x^) ^^"^ arccosh(x/?/), x > y, 
and E{x) is the complete elliptic integral of the second kind.'^Sl 
The constant coefficient £o(C) can be written as the sum 

50(0=^^(0+4 



of a direct ("ring-diagram") term 6^(0 and an exchange term Eq. Following Onsager's worlP^ 
on the three-dimensional gas, the exchange term was found by Isihara and loriattP^ to be 



(15) 



4 = /3(2)--^/3(4) = +0.114357, 

TT 



(16) 
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FIG. 1. T_2(C), Tr-i(C), ^^0(0, Tfj(C) and Tt,(C,) as functions of C- 

where /5 is the Dirichlet beta functioiP^ and G = /3(2) is Catalan's constant. We note that 
£q is independent of C, and the spin-scahng function therefore takes the trivial form 



(17) 



The direct term has not been found in closed form, but we now show how this can be 
achieved. Following Rajagopal and Kimball/^ we write the direct term as the double integral 



^o(C) 



OO /"OO 



oo -'0 



Q 



U 



dq du, 



where 



(18) 



lu — i] m + 1 

2 A2 



(19) 



In the paramagnetic (C = 0) case, the transformation s = (f /A — v? and t = qu yields 



em) 



1 

'2^ 



OO fOO 



OO JO 



+ t2 



IS - 1)2 + t2 + s _ 1 
Vs2 + t2 + s 



1/2" 



n 2 



rftrfs, (20) 
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and, if we adopt polar coordinates, this becomes 



1 



OG PIT 




1 

2^ 



\/l — 2r cos 9 + r"^ — 1 + r cos 9 
r(l + cos 6*) 

9 



d9dr 



Q 

2 In 2 - (tt - 6*) tan - - 2tan^ - In ( sin 



d9 



= ln2 - 1 

= -0.306853..., 

which confirms Seidl's numerical estimat^^ 



£^(0) = -0.30682 ±0.00012. 



In the ferromagnetic (C = 1) case, Eq. (18) yields 



a / N 1 » / N ln2 — 1 
4(1) = 2^(0) = = -0.153426 .... 

In intermediate cases, where < C < 1, we define the spin-scaling function 

4(C) 



T^(C) 



and, from (18), we have 



T^(C) 



1 



1 



2 47r(ln2-l)yo 



em' 



°° '"^ idz 

Pk^{r, z)Pk^{r, z) dr, 




where 



^, , ^ \/rz -k"^ + Jrlz-k^ 
Pk{r, z) = l ^ . 

v^(v^+ 1/v^) 



Integrating over r gives 



T^(C) = 

where 

Lkf,kj, (z) = -k^ In k^ - ki In % 
1 



1 



1 



2 47r(ln2 - 1) 



^ .Adz 

^k^,k^[z) , 



+ 



[z + iy 



{zk-i - ki) \n{zk^ - k^) + {zk^ - k^) \n{zki - k^) 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



it: 



{kf - 2zk-^k^ + kf) + 2z{k^ + k^f \n{k^ + k^) - z{zk^ - 2k-^ki + zkf) In , 



(28) 
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TABLE I. Energy coefficients and spin-scaling functions for 2-jeIlium in tlie high-density limit. 



Term 



Coefficient 



e{0) 



e{l) 



T(C) 



r, Inr, 



^-2(0 

^-i(C) 
eo(C) 



1 

2 

3^ 
ln2 - 1 

m - 



3tt 



37r 
ln2 - 1 
2 

m - 

1/10 



4 V 37r 



Eq. ((5j) 
Eq. (|6j) 



Eq. (|29|) 
1 



Eq. (12) 



and contour integration over z eventually yields 



1-C 



2 4(ln2-l^ 



21n2 - 1 



1 + C 
1-C 



(29) 



This is plotted in Fig. [l] and agrees well with Seidl's approximation,'^ deviating by a 
maximum of 0.0005 near ( = 0.9815. 

In conclusion, we have shown that the energy of the high-density spin-polarized two- 
dimensional uniform electron gas can be found in closed form up to 0{rs). We believe that 
these new results, which are summarized in Table |H will be useful in the future development 
of exchange-correlation functionals within density-functional theory. 
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